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£f) • In this paper, we consider some frustrated spin models for which the ground states are 

known exactly. The concurrence, a measure of the amount of entanglement can be calculated 
exactly for entangled spin pairs. Quantum phase transitions involving macroscopic magneti- 
zation changes at critical values of the magnetic field are accompanied by macroscopic jumps 
in the T = entanglement. A specific example is given in which magnetization plateaus give 
rise to a plateau structure in the amount of entanglement associated with nearest-neighbour 



(N 
> 



O 

00 



s 

T3 



Abstract 



bonds. We further show that macroscopic entanglement changes can occur in quantum phase 
transitions brought about by the tuning of exchange interaction strengths. 
PACS Nos.: 03.65.Ud, 75.10.Jm 



Entanglement is a characteristic feature of quantum mechanical systems which has no classical 
analogue [1]]. The state of a pair (or more than a pair) of quantum systems is entangled if the 
corresponding wave function does not factorize, i.e., is not a product of the wavefunctions of the 
individual systems. A well-known example of an entangled state is the singlet state of two spin— | 
particles, ^= (|tl) — lit))? which cannot be written as a product of the spin states of individual 
spins. Measurement on one component of an entangled pair fixes the state of the other implying 
non-local correlations. Interest in quantum entanglement is extensive because of its fundamental 
role in quantum communication and information processing such as quantum teleportation (2|, 
superdense coding ||, quantum cryptographic key distribution [|J etc. Experimental implemen- 
tation of some of the protocols has so far been achieved in simple physical settings. Solid state 
devices, specially, spin systems have been proposed as possible candidates for large scale realiza- 
tions H |||. In particular, the Heisenberg spin-spin exchange interaction gives rise to entangled 
states in spin systems and has been shown to provide the basis for universal quantum computation 
[|7|, U . Examples of other interacting many body systems in which entanglement properties have 



been studied include the harmonic chain ||, the ID Kondo Necklace model QT0|| and the BCS 
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Entanglement in a state like its energy is quantifiable and has been computed both at T = and 
at finite T (thermal entanglement) for a variety of spin models in both zero and finite external 
magnetic fields. The models include the Heisenberg XX, XY, XXX, XXZ and transverse 
Ising models in one dimension (ID) [12-18]. The computational studies show that the amount of 
entanglement between two spins in a multi-spin state can be modified by changing the temperature 
and/or the external magnetic field. Some recent studies have explored the relations between 
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entanglement and quantum phase transitions in the ferromagnetic (FM) XY model in a transverse 



magnetic field and in a special case of the model, namely, the transverse Ising model |16|, [17 
A quantum phase transition (QPT) can take place at T = by changing some parameter of 
the system or an external variable like the magnetic field |l9fl. The ground state wavefunction 
undergoes qualitative changes in a QPT and model studies indicate that entanglement develops 
special features in the vicinity of a quantum critical point. 

The studies carried out so far have been confined to ID systems with only nearest-neighbour 
(NN) exchange interactions. There are several quasi-lD and 2D antiferromagnetic (AFM) spin 



models which describe frustrated spin systems |g0|. Frustration arises if conflicting interactions 



are present in the system, i.e., when all the interactions between spins cannot be simultaneously 
satisfied. A good example of a frustrated spin system is the AFM Ising model on the triangular 
lattice. An elementary plaquette of the lattice is a triangle. The Ising spin variables have two 
possible values, ±1, corresponding to the up and down spin orientations. An antiparallel spin pair 
has the lowest interaction energy — J. A parallel spin pair has the energy +J. In an elementary 
triangular plaquette, there are three interacting spin pairs. Due to the topology of the plaquette 
(odd number of NN bonds), all the three spin pairs cannot be simultaneously antiparallel in the 
ground state. There is bound to be one parallel spin pair giving rise to frustration in the system. 
Consider another example in which the spins in a linear chain interact via NN as well as next- 
nearest-neighbour (NNN) AFM interactions. In a triplet of successive spins, the two NN spin pairs 
and one NNN spin pair cannot be simultaneously antiparallel and the linear spin chain is frustrated. 
Frustration in a system may occur due to the topology of the lattice (triangular, kagome etc.) or 
due to the inclusion of further-neighbour interactions leading in most cases to spin-disordered 
ground states. The models exhibit QPTs as the exchange interaction parameters are tuned to 
particular ratios or at critical values of the external magnetic field. Some of the models exhibit 
the phenomenon of magnetization plateaus in which plateaus appear in the magnetization/site m 
versus the external magnetic field h curve at quantized values of m (m and h are chosen to be 
dimensionless) |21||. The condition for the appearance of plateaus is 



Su — mu = integer (1) 

where Su and mu are the total spin and magnetization in unit period of the ground state. The 
plateaus indicate that the spin excitation spectrum is gapped so that the magnetization remains 
constant. In this paper, we consider some spin models for which the ground states are known 
exactly in certain parameter regimes. An exact measure of the entanglement between two spins 
can be obtained in these states in a straightforward manner because of the simple structure of 
the ground states. We report two major results. In the presence of an external magnetic field, 
macroscopic magnetization jumps at critical values of the magnetic fields (first order QPTs) are 
accompanied by macroscopic changes in the amount of pairwise entanglement. Furthermore, some 
examples are given in which the entanglement structure is modified due to QPTs (again first order) 
brought about by the tuning of exchange interaction strengths. 

A measure of entanglement between the spins A and B is given by a quantity called con- 



currence [12], [0|. To calculate this, a knowledge of the reduced density matrix pab is re- 
quired. This is obtained from the ground state wavefunction by tracing out all the spin de- 
grees of freedom except those of the spins A and B. Let pab be defined as a matrix in the 
standard basis {|TT) > IT I) ; I IT) > 111)}- O ne can define the spin-reversed density matrix as p = 
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{tjy <£> cr y ) p* (a y <8> cr^), where cr y is the Pauli matrix. The concurrence C is given by C = max{Ai — 
A2 — A3 — A 4 , 0} where X^s are the square roots of the eigenvalues of the matrix pp in descending 
order. The spins A and B are entangled if C is nonzero, C = implies an unentangled state 
and C — 1 corresponds to maximum entanglement. The models we consider in this paper belong 
to the Majumdar-Ghosh (MG) |22| and the Affleck-Kennedy-Lieb-Tasaki (AKLT) families of 
models. The MG model is the simplest frustrated model in ID. The spins have magnitude \ and 
the Hamiltonian is given by 

N j N 

Hmg = J 2^ ^i- S i+ i + — Si. S i + 2. (2) 
i=i i=i 

N is the total number of spins and the strength of the next-nearest-neighbour (NNN) exchange 
interaction is half that of the NN interaction. The boundary condition is periodic. The exact 
ground state is doubly degenerate with the wave functions 



[12] [34] [N — IN] 

[23] [45] [Nl] 



(3) 



where [Im] denotes a singlet (valence bond (VB)) of two spins located at the lattice sites I and 
m. Thus in the ground state NN spin pairs in singlet spin configurations are maximally entangled 
with concurrence C — 1. The value of C for all the other NN or further-neighbour pairs is zero. 
This is consistent with the special property of entanglement that in a set of three spins A, B 



and C, if A and B are maximally entangled, the entanglement between A and C is zero (12 
Translational invariance in the ground states can be restored by taking linear combinations of the 
ground states, < ^ 1 ^ 2 . In this case, all NN pairs are entangled with C = 0.5. The value of C in the 
case of the S = ~ Heisenberg AFM spin chain is C = 0.386 
the NN entanglement in a spin system. 



12 . Frustration appears to increase 
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FIG. 1: A two-leg frustrated spin ladder. The spin-spin exchange interaction strength along the 
rung is Jr. The inter-leg NN and the diagonal exchange interactions are of equal strength J. 

We now give examples of some spin models in the presence of an external magnetic field for 
which macroscopic magnetization jumps are accompanied by macroscopic jumps in entanglement. 
The first model is a frustrated two-leg S - 
is given by 



ladder model |24j] (Fig. 1) the Hamiltonian of which 



N 



(ij) 



Si.Sj-hJ^St 



(4) 
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The exchange interaction strength Jij = Jr along the rungs. The intra-chain NN and the diagonal 
exchange interactions are of equal strength J. When h = and > A (A ~ 1.401), the exact 
ground state consists of singlets (VBs) along the rungs. The spin pairs along the rungs are thus 
maximally entangled and all other spin pairs are unentangled. We define a quantity 



(5) 



which is the average concurrence per rung of the ladder. The summation is over the rung index 



with y being the total number of rungs. Since all the other spin pairs are unentangled, the average 
concurrence per NN bond, fwN = § (the NN bonds include y rung bonds and N intra-chain NN 
bonds). The magnetization properties of the frustrated ladder model are simple [25|]. Fig. 2(a) 
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fa) (b) 
FIG. 2: The frustrated spin-ladder model (Fig. 1) exhibits plateaus in (a) magnetization/rung M 
versus an external magnetic field h and (b) average concurrence per rung / versus h. 



shows the magnetization per rung M as a function of the magnetic field h. For < h < h ci — Jr, 
M = 0. For h Cl < h < h C2 = Jr + 2J, the rungs are alternately in singlet and S z = 1 triplet spin 
configurations in the ground state. The value of M is now |. For h > h C2 , saturation magnetization 
is obtained with M — 1. One can verify that the quantization condition (1) is obeyed at each 
plateau. Fig. 2(b) shows the plot of the average concurrence per rung / versus the magnetic field 
h. A similar plot is obtained for f^^ with the difference that /nn = | for < h < h Cl and 
Inn = | for h Cl < h < h C2 . For h Cl < h < h C2 , half the total number of rungs are in S z = 1 
triplet spin configurations in which the two spins of a rung point up. The parallel spin pairs are 
unentangled with C = 0. The other rungs are in singlet spin configurations so that f =\- In the 
fully polarised state / = 0. Fig. 2(b) provides an example of a changing magnetic field giving rise 
to macroscopic changes in entanglement. The ground state of the frustrated spin ladder model 
for h = has a simple structure and can be expressed as a product over rung singlet states. The 
reduced density matrix pab, needed to calculate the entanglement between the spins A and B, 
can be calculated exactly and analytically in a straightforward manner. Once the matrix elements 
of pab are known, the concurrence C can be calculated. Because of the uncomplicated structure 
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of the ground state, C can alternatively be calculated using simple arguments. Each spin in the 
ladder belongs to a rung singlet or VB so that it is maximally entangled with its partner spin. 
Thus, the entanglement between a spin j and any other spin, belonging to other rung singlets, is 
zero. The concurrence C of a maximally entangled pair is 1 so that /, the average concurrence 
per rung, is 1 for < h < h Cl . Similar arguments hold true for h ci 
a different ground state, again with a simple structure. 



< h < h C2 when the ladder has 
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FIG. 3: Two-chain ladder model consisting of four-spin plaquettes coupled to two-spin rungs. 

The exchange 
interaction strengths are as shown in the Figure. 



Fig. 3 shows another spin ladder model with modulated exchange interactions |f2^]. The 
model consists of four-spin plaquettes coupled to two-spin rungs (solid lines) through NN and 
diagonal exchange interactions (dotted lines) of strength J 2 . Within each plaquette, the NN 
and diagonal exchange interactions are of strength J% and J3 respectively. The rung exchange 
interaction strength is Jr. In a wide parameter regime, the exact ground state has a simple product 
form: the plaquettes and the rungs are individually in their ground state spin configurations. The 
ground state of a plaquette is a resonating valence bond (RVB) state with wavefunction given by 
iPrvbi{iPrvb2) for J3 < Ji{J% > Ji). The ground state of a rung is a singlet. The RVB states are 
given by 



!'RVBl(RVB2) — 



+ 

(-) 



(6) 



The solid lines represent singlets and the arrow signs are drawn according to the phase convention 
that in a VB between the sites i and j, if the arrow points away from the site i, then the spin 
configuration is ^= (|f (i) j. (j)) — || (i) f (j))). It is easy to check that the NN spins in \iPrvbi) 
are entangled with concurrence C = 0.5. The NNN spins along the diagonals are unentangled. On 
the other hand, the NN spins in \1pRVB2) are unentangled and the NNN spin pairs are entangled 
with C = 1. At Ji = J 3 , there is a QPT from the ground state in which the plaquette spin 
configurations are described by \4>rvbi) to the ground state in which the same are described by 
\ipRVB2)- In both the phases, the rungs are in singlet spin configurations. The QPT is accompanied 
by macroscopic changes in the amount of NN and NNN entanglements. The average concurrence 
per NN bond, /nn, in the full ladder is ~ for J 3 < J\ and | for J 3 > J\. The average concurrence 
per NNN bond, fo, is for J3 < J\ and | for J3 > J\. In a finite magnetic field h ^ 0, the 
exact ground state maintains its product form in an extended parameter regime. This gives rise 
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to magnetization plateaus in the magnetization/site m versus h curves. Again, the jumps in the 
magnetization are accompanied by jumps in the amount of entanglement. To give one specific 
example, consider the case J3 < J\. The average concurrence/NNN bond in the full ladder is zero 
for < h < h ci — J\ and \ for h Cl < h < h C2 = 2J 2 . At h Cl , m jumps from zero to the value |. 
Each plaquette is in the ground state spin configuration | (|Tttl) — ITTIT) — ITITT) + IITTT)) m 



the field range h ci < h < h C2 . Schulenberg et al. [27| have constructed exact eigenstates which 




FIG. 4: Kagome lattice with N spins. The circles mark the g- hexagons in which independent 

magnons can be localised. 

consist of independent, localized one-magnon states in a class of frustrated spin lattices and have 
shown that these are the ground states in high magnetic fields. The eigenstates are obtained 
provided certain conditions are satisfied. In a kagome lattice, the magnon states are localised in 
Y hexagons where N is the total number of spins. The hexagons in which the magnon excitations 
occur are isolated from each other (Fig. 4). Above the saturation magnetic field, all the spins in 
the lattice are in a FM spin configuration |0) and each spin pair is unentangled. The state |0) is 
the vacuum for magnon excitations. At a critical field below the saturation magnetic field, the 
exact eigenstate consisting of ^ localised one-magnon excitations becomes the ground state. The 
localised magnon excitation has the wave function 

|i) = ^E(-i)^rio). (7) 



1=1 



This gives rise to a macroscopic change in the magnetization curve. In each of the ^ hexagons, 
a spin is equally entangled with all the other five spins and the magnitude of the concurrence is 
C = |. This signifies a macroscopic change in the amount of entanglement. The magnetization 
jump occurring at a critical value of the magnetic field signifies a first order QPT. We have 
shown through specific examples that such QPTs give rise to macroscopic jumps in the amount of 
entanglement associated with NN and/or further-neighbour spin pairs. 

QPTs can also be brought about by tunning exchange interaction strengths. We have already 
given an example of this in the case of the ladder model shown in Fig. 3. For h — 0, a QPT 
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occurs at 



Ji 



1. In the frustrated two-leg ladder model shown in Fig. 1, a QPT takes place at the 
critical value {^f) c — 1.401. Below the critical point, the ground state is that of an effective S = 1 
chain with the 5* = 1 spins forming out of the pairs of S — \ rung spins. The spin ladder is now 
in the Haldane phase of a S — 1 chain. Above the critical point, the ladder is in the rung-singlet 



(RS) phase in which the rung spins are in singlet spin configurations. Kolezhuk and Mikeska [28 



have proposed a generalised frustrated ladder model in which the first order QPT between the 
RS and the Haldane phases can be studied in an exact manner. This model includes biquadratic 
interactions besides NN and NNN (diagonal) ones. The Hamiltonian is given by, 



H — hi t i + i, 



yi. 



hi,i+\ — — ( S 1,1. S 2,i + S S 2,i+l) + ( S l,i- S 1,1+1 + S 2,i- S 2,i+l) 

J rV2\ S S 2,i+l + S 2,i- S l,i+l) + X\{ S i,j. S l,i+l)( S 2,i- S 2,i+l) 



+X2( S i,j. S 2,i+l)( S 2,i- S l,i+l) 



(8) 



4 4 
where, x x = -(3 - 2y 2 ), a;2 = 7(3y 2 - 2) 



(9) 



The indices 1 and 2 correspond to the lower and upper legs of the ladder respectively and i is 
the rung index. The exact phase boundary between the RS and the Haldane phases is given by 
Hi = |(1 + y 2 ). The ground state in each phase is known exactly. In the Haldane phase the ground 
state is the AKLT state p3|| . The ground state energy/rung in the RS and the Haldane phases 
are E RS = -jUi + ^(1 + y 2 ) and E AKLT = \y x - + y 2 ) respectively. 

As mentioned before, in the RS phase the spin pairs along the rungs are perfectly entangled 
with concurrence C = 1. All the other spin pairs are unentangled. The amount of entanglement 
in the AKLT phase can be computed in the following manner. For [H, S z ] = 0, where S z is the z 
component of the total spin, the reduced density matrix of a spin pair located at the sites % and j 
has the form |T^] 

/ u + 
wi z 
Z W2 




Pij 



\ 



\ 



u. 



(10) 



/ 



The concurrence quantifying the entanglement is given by 



C = 2max[0, 



(11) 



Wang and Zanardi |^] have shown that the matrix element of can be expressed in terms of the 
various correlation functions G Q/ g = ((Ji a <Jjp) = Tr{a ia ajpp), (a = x,y,z), where p is the density 
operator. The magnetization/site m = jfTr(^2 i=1 crfp). In particular, the following relations hold 
true: 

1 

G zz ) 



u± 



-(1 ±2m 
4 V 



7 



(G xx + G 



mi I 



(12) 



The correlation functions in the AKLT phase can be calculated in an exact manner using the 
transfer matrix method in the matrix product (MP) formalism |3IJ. The AKLT ground state can 
be written in the MP form as 



N_ 

2 



AKLT 



TrJl 



9i 



with, gi 



l*b>, 



-V2|f+>, 
- \to) t 



(13) 



(14) 



The product in Eq. ( |i~3"D is over the y rungs of the ladder. The two spin-|'s of each rung are in 
a triplet spin configuration in the AKLT phase giving rise to an effective spin 1. In Eq. flUl) , the 
states \tfj) with /x = +1, and —1 represent a spin state with S z = +1, and —1 respectively. 
Calculation of the correlation functions G aa (a = x, y, z) of the spin ~ pairs in the AKLT state 
can be carried out following standard procedure 
correlation functions are: 



30 . We quote the final results. The various 



\Cl,na cr l,n+la) — (&2,na&2,n+la) — \0'l,na°'2,n+la) — (01,n+la°2,na) 



= -q (a = (15) 

(ffl.na^.no) = g (a = x,y,z). (16) 

The correlation functions, in which the distance I separating the rungs on which the spins are 
located is > 1 (in Eq. ( [Of ) 1=1), involve the factor 4(— 1) J 3~' -1 . With the knowledge of the 
correlation functions, the concurrence can be determined using the relations ( |ll~l ) and ([12] ) (m in 
the AKLT phase is zero). One finds that in the AKLT phase, the spin pairs along the rungs are 
unentangled. In the RS phase, the same spins are maximally entangled. The intra-leg NN and 
the NNN (diagonal) spin pairs are entangled in the AKLT phase with concurrence C = | in each 
case. In the RS phase, these spin pairs are unentangled. Again, the first order QPT from the RS 
phase to the AKLT phase is accompanied by macroscopic changes in the entanglement structure. 
The generalised ladder model studied by Kolezhuk and Mikeska has a rich phase diagram 



28| describing both first order and second order QPTs. At the second order phase boundary, 
the gap in the excitation spectrum goes to zero. There are five possible phases: RS, AKLT, 
FM, Dl and D2. The phases Dl and D2 have spontaneously broken symmetry but the exact 
ground states are not known in these phases. The phase boundaries separating the FM-D1, 
AKLT-D1, AKLT-D2, RD-D1 and RD-D2 are known exactly and have been determined in the 
MP formalism. The corresponding QPTs are second order transitions. The two first order phase 
boundaries separating the RD-AKLT and AKLT-FM phase are also known exactly. Preliminary 
calculations |3T] near the second order phase boundaries suggest that the pairwise entanglement 
does not extend beyond the NNN distance. The same is true in the case of first order QPTs. 
Osterloh et al. [TH| have considered the range of entanglement in the vicinity of the quantum 



8 



critical point of the transverse Ising model in ID. They have found that even at the critical point, 
where the spin-spin correlations are long ranged, the concurrence is zero for spin pairs separated 
by more than NNN distance. There is, however, one crucial difference between this model and 
the AKLT-type models described by MP states. The MP states are finitely correlated, i.e., the 
spin-spin correlation function (S£S* +l ) = A s e « decays exponentially with the correlation length 
£ equal to a few lattice spacings. As the transition point r = t c (r is some model parameter) 
is approached, the correlation length £ either does not exhibit any singularity or diverges in a 
power-law fashion as (r — r c ) 2 . In the latter case, however, the prefactor A s oc (r — r c ) becomes 
zero at the transition point [28]. Thus long range spin-spin correlations cannot develop in the 



system. 

In summary, we have shown through specific examples that first order QPTs can bring about 
macroscopic changes in the amount of pairwise entanglement in spin systems. We have given the 
specific example of a spin model in which magnetization plateaus give rise to a plateau structure 
in the average concurrence per rung as well as NN bond. There are several low-dimensional 
AFM compounds in which magnetization plateaus have been observed experimentally [[32|]. The 
appearance of magnetization plateaus has been explained in terms of metal-insulator transitions 
of magnetic excitations driven by a magnetic field p3| . In the insulating (plateau) phase, the 
magnetic excitations give rise to crystalline order and in the metallic (non-plateau) phase they 
are itinerant. It will be of interest to explore the possibility of a plateau structure in the amount 
of T = as well as finite T entanglements in such systems. Our study shows that an external 
magnetic field can be employed to give rise to large changes in the amount of entanglement and 
provides the basics for the construction of an entanglement 'amplifier' or 'switch'. There is a large 
number of spin models which exhibit QPTs of significant interest. Some of these models describe 
2D systems. Examples include the Shastry-Sutherland model fl33J, the S — % AFM model on the 
|-depleted square lattice [[35|], a lattice of weakly-coupled two-leg ladders |36| etc. Some of these 



models are reviewed in Refs. pO], |37|| . All these models exhibit second order QPTs from a spin- 
disordered gapped phase to a gapless phase with long range spin-spin correlations. Osterloh et al. 



16| and Osborne and Nielsen [T7|] have found evidence of entanglement showing scaling behaviour 
in the vicinity of the quantum critical point. Similar studies should be carried out for the spin 
models mentioned to augment our knowledge of the relationship between QPTs and entanglement. 
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